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Prescribed  by  ANSI  Std  Z39-18 


Fuel  Economy,  Mobility  and 

Safety 


System  Level:  Battery,  Gearbox,  Occupant  Compartment  Design;  Motor  Map  Selection 


Selected  Motor  Map,  System  Level 


Resizable  Motor  Map,  Subsystem  Level 
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Motivation 


Underbody  blast  events  are  a  top  threat 
facing  U.S.  Army  ground  personnel 


?O07  ^  '?GOp.  3% 


LED  Deaths  By  Year 

!■ 

Hostile  Deaths 
■  1ED 


300 


200 


ISO 


http://www.focusblog.ro/wp-content/uploads/2010/03/LAND_M1114_HMMWV_IEDed_lg.jpg,  accessed  April  29,  2010 
jCasuahies(2()10).“IED  Fatalities.”  http://icasualties.org/oef,  accessed  April  6^2010. 
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Motivation 


Vehicle  weight  has  mixed  effects  on 
different  design  objectives 


High  Mobility  Multipurpose  Wheeled 
Vehicle  (HMMWV) 

2,700  kg 

http://www.amgeneral.com/vehicles/hmmwv/a2-series/details/ml097a2-base 

http://www.globalsecurity.org/military/systems/ground/caiman-specs.htm 


Mine  Resistant  Ambush  Protected 
Vehicle  (MRAP) 

14,000  kg 
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Research  Objective 


Multi- objective  optimization  of  ground 
vehicles  for  reduced  weight  and 

occupant  injury 

Determine  occupant  injury  as  a 
response  to  structural  and  occupant 
compartment  design  parameters 


Develop  surrogate  models  for  vehicle 
and  occupant  responses  to  a  blast  event 


Account  for  uncertainty  in  blast 
location  and  size 


http://c0378172.cdn.cloudfiles.rackspacecloud.com/7770_9080764544.jpg,  accessed  on  April  29,  2010. 
http://mocoloco.com/art/archives/pickering  land  mine  mar  06.jpg,  accessed  on  April  14,  2010. 
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Optimal  Design  Laboratory 


Modeling  Approach 


I  nputs: 

Vehicle  Mass 

Charge  Location  (x,  y  coordinates) 
Charge  Mass 


Underbody  Blast  Simulation 


Outputs: 

Upper  Neck  Axial  Force 
Lower  Lumbar  Axial  Force 
Lower  Tibia  Axial  Force 


Drop  Tower  Simulation 


Blast  Pulse  of  Vehicle 


4 


I  nputs: 

Blast  Pulse  (magnitude  &  duration) 
Seat  Cushion  Stiffness 
Seat  Energy-Absorbing  (EA) 

System  Stiffness 
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Optimal  Design  Laboratory 


Charge  Uncertainty 


Field  data  about  charge  distribution  is 
sensitive,  so  I  postulate  distributions: 


Charge  longitudinal/ 
x-location  ~  U(a,b) 
(m) 


Charge  lateral/ 
y-location  ~  U(a,b) 
(m) 


Charge  mass  ~  N(fj,cr) 
(TNT-equivalent) 
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Uncertainty  in 
Vehicle  Peak 

Acceleration 

_ 

Itb^ . 

Charge  Uncertainty 
Propagation 


Drop  Tower  Simulation 


Uncertainty 
in  Occupant 
Body  Forces 


Upper  Neck 
Compression 


4 


Lumbar  Lower  Leg 
Spine  Compression 


Compression 

Uncertainty  in 
Optimal  Seat  System 
Design  Variables 


Optimization  to 
Minimize  Injury 
Probability 

4 
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Structural  Model 


I  nput  Variables: 

Vehicle  Mass  (mv) 
Charge  Location  (xc,  yc) 
Charge  Mass  (mc) 


Output: 

Blast  pulse  (aDeak) 


Time  (ms) 


Surrogate  model  from  linear  regression  on  100 

data  points: 

1  xc  yc  mc 

apeak  =  Po  +  Pi  ~  +  ihxc  +  PiYc  +  p4mc  +  ih  ~  +  ft>  “  +  Pi  ~  +  Ps)7cmc  +  p9>rc 


V 


Livermore  Software  Technology  Corporation  (2007).  LS-DYNA  Keyword  User’s  Manual,  http://lstc.com/pdf/ls-dyna 971 manual k.pdf,  accessed  April  29,  2010. 
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-500 


Blast  Pulse  Uncertainty 


Peak  Acceleration  (g's) 


Peak  accelerations  for 
4,000  kg  vehicle 


Mean 


0  2000  4000  6000  SOOO  1QOOO  12000  14000 
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Distribution  moments  plotted  versus  vehicle  mass 
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Occupant  Model 


I  nputs: 

Blast  Pulse  (apeak) 

Seat  Cushion  Foam  Stiffness  (sc) 
Seat  EA  System  Stiffness  (sEA) 


Outputs: 

Upper  Neck  Axial  Force  (Fneck) 
Lower  Lumbar  Axial  Force  (F|Umbar) 
Lower  Tibia  Axial  Force  (Ftibia) 


Arepally,  S.  et.  al.  (2008).  Application  of  Mathematical  Modeling  in  Potentially  Survivable  Blast  Threats  in  Military  Vehicles. 
http://www.dtic.mil/cgi-bin/GetTRDoc?AD=ADA496843&Location=U2&doc=GetTRDoc.pdf,  accessed  on  April  29,  2010. 
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Occupant  Model 

Surrogate  model  from  linear  regression  on  500 

data  points: 

p  _  f  3.809— 0.03934s  EA +0.4289sc  +  0.003446fip  eak  +0.0002 16'ldfg^  dpeak~ 0.000001781rtpea^^ 

*  neck  ~ 

F spine  =  383  -  462 sEA  +  416sc  +  1.4 apeak  +  262 sEAsc  +  0.7 sEAapeak  +  scapeak  -  232 sc2  -  Q.OOGba^ 
^combined  tibia  07  +  63S£^  495sc  +  3.7 Qpeak  0.1 6sEACLpeak  0.38scCJpefl^  +  99sc  +  0.0003tJpen^. 

U.S.  Army  aims  for  no  more  than 
10%  probability  of  moderate  injury 
(AIS2+) 

Research  and  Technology  Organisation  (2007).  “Test  Methodology  for  Protection  of  Vehicle  Occupants  against  Anti- Vehicular  Landmine 
Effects.”  North  Atlantic  Treaty  Organisation,  Neuilly-sur-Seine  Cedex,  France.  Accession  number  RTO-TR-HFM-090. 


Thresholds: 

$ 

^neck  4  kN 

^lumbar  —  6. 7  kN 

<1 

Ftibia  =  5.4  kN 
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Optimization  Formulation 

General  Safety  Objective:  minimize  occupant  injury 


What  is  the  explicit  objective  function? 


minf(x)  =  probability  of  AIS2  Injury 

Complication:  unknown  injury  probability 
distributions 


Thresholds: 

ii  - _ 

^neck  4  kN 

f) 

I 

^lumbar  —  6. 7  kN 

Ftibia  =  5.4  kN 

min  f(x)  =  body  forces  experienced  when  vehicle  is  attacked 

Complications:  uncertainty  in  charge  parameters, 
multiple  body  forces  of  interest 

-500  0  500  1000  1500  2000  2500 

Peak  Acceleration  (g's) 
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Formulation  1:  Model 

Objective:  minimize  the  maximum  of  the 
body  forces  (percentage  of  threshold) 


minSM,Sc  max  | 

( Fneck  F lumbar  Ftibia(combined)\ 

{  4 kN  '  6.7 kN  '  5.4 kN  ) 

where 

a peak  l^OLpeajcC^'v')  —  4  X  1 0  I7lv 

Fneck  ^neck\sEAr  scf^peak) 

F 'lumbar  ^lumbar \sEAf  sct  ^peak) 

F combined  tibia  ^ combined  tibia  EA>  sc 

subject  to 


lb  <  sEAfsc  <  ub 


UNCLASSIFIED 


14 


Optimal  Design  Laboratory 


- Direction  of  Improved  Occupant  Safety 

Maximum  Percentage  of  Injury  Threshold 


Formulation  1:  Results 


Objective:  minimize  the  maximum  of  the 
body  forces  (percentage  of  threshold) 


0  2000  4000  6000  8000  10000  12000  14000 

Total  Vehicle  Mass  (kg) 

< -  Direction  of  Improved  Mobility  and  Range 


UNCLASSIFIED 


15 


Formulation  1  Limitation 

Minimizes  body  forces  for  a  given  vehicle 
mass  for  50th  percentile  of  charges 


Peak  Acceleration  (g's) 
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Optimal  Design  Laboratory 


Formulation  2:  Model 

Objective:  minimize  the  probability  of 
“failure"  to  meet  injury  threshold 

Pf  ~  1  —  (Apeak) 

!  _|_  erf  |  “PeaX-^peak 
\  J 2(Tapeak 

=  4  x  106m-1023 

%ea^mu)  =  2  x  106m“1035 

subject  to  Fneck  =  g1(sEA, sc,  apeak)  <  4000 

Peak  Acceleration  (g's) 

Plumb ar  9z\^EA>  ^ct^-peak)  —  0700 

Pcombined  tibia  93\SEA>  sc>  apeak J  —  o400 

lb  <  sEA,sc  <  ub 


where 


^  (Apeak)  7 
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Probability  of  Failure 


Formulation  2:  Results 


Objective:  minimize  the  probability  of 
“failure"  to  meet  injury  threshold 


Failure  Probability  vs.  Vehicle  Mass 


0  2000  4000  6000  8000  10000  12000 

Vehicle  Mass  (kg) 


Vehicle 
Mass  (kg) 

Probability  of 
Failure 

2000 

4.60E-01 

2500 

2.45E-01 

3000 

9.93E-02 

3500 

2.97E-02 

4000 

6.43E-03 

4500 

9.90E-04 

5000 

1 .07E-04 

5500 

8.06E-06 

6000 

4.20E-07 

6500 

1.51E-08 

7000 

3.69E-10 

7500 

6.16E-12 

8000 

6.99E-14 

8500 

5.55E-16 

9000 

0.00E+00 

SEA  ~  l-^c  ~  2.0 
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Occupant  Model  with  Floor  Pad 


Inputs: 

Blast  Pulse  (apeak) 

Seat  Cushion  Foam  Stiffness  (sc) 
Seat  EA  System  Stiffness  (sEA) 
Floor  Pad  Foam  Stiffness  (sf) 


Surrogate  model  from  linear 
regression  on  300  data  points: 


U  —  A.  3-84  +  0.12sEj4  +0.88sc  +  0.002ap^afc +0. 0585^  Sc  +  0. 000084%^  apeajt-0.000063scapgajt-0.058s|^-0.14s|-0.00000054a^eafc] 
**neck  ~ 


F  ■ 

1  spine 


(s.664  +  0.12s£^  +  0.81Sc  +  0.002apeafc+0.062sE>1  SC  +  0.QQQQ87 SEAapBak-0.QQQQ6Qscapeak  -0.039s^A -0.13s| -0.000000 S6apeafc) 


^combined tibia  =  332  -  245sc  -  80.23sr  +  1.3 apeak  +  35.84 scsf  +  14.0 sf  +  0.0012 a^eak 


Optimal  Design  Laboratory 


UNCLASSIFIED 


19 


- Direction  of  Improved  Gccupa nt Safety 

Maximum  Percentage  of  Injury  Threshold 


Results  with  Floor  Pad 


Objective  1:  minimize  the 
maximum  of  the  body  force 
percentages 


Objective  2:  minimize  the 
probability  of  "failure"  to 
meet  injury  threshold 


Occupant  Safety  versus  Vehicle  Mass 


0  2000  4000  6000  S00Q  10000  12000  14000 

Total  Vehicle  Mass  (kg) 

^ -  Direction  of  improved  Mobility  a  rid  Range 


Failure  Probability  vs.  Vehicle  Mass 


0  2000  4000  6000  S000  10000  12000 

Vehicle  Mass  (kg) 
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Summary 

■  Developed  a  modeling  approach  to  evaluate 
structural  and  occupant  responses  to  ground 
vehicle  underbody  blasts 

■  Fit  surrogate  models  to  reduce  computational 
expense 

■  Demonstrated  two  optimization  formulations  and 
their  results 

D  Accounted  for  uncertainty  in  charge  parameters 

D  Quantified  negative  correlation  between  vehicle  mass 
and  occupant  injury  probability 

■  Added  floor  padding  to  reduce  tibia  impact 
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Ongoing  Work 


tank-net.org,  accessed  Apri  27,  2010.  ROllOVer  SafetV  mOdelinO 

http://www.usaasc.info/alt  online/images/080901  Photo2.jpg,  accessed  April  29, 2010.  J  u 
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I^H||  Optimal  Design  Laboratory 


Motivation 


Underbody  blast  events  are  a  top  threat 
facing  U.S.  Army  ground  personnel 


IED  Fatalities 


Period 

IED 

Total 

Pet 

2001 

0 

4 

0.00 

2002 

4 

25 

16.00 

2003 

3 

26 

11.54 

2004 

12 

27 

44.44 

2005 

20 

73 

27.40 

2006 

41 

130 

31.54 

2007 

78 

184 

42.39 

2008 

152 

263 

57.79 

2009 

275 

450 

61.11 

2010 

86 

150 

57.33 

UNCL4§^aWl0X 


“IED  Fatalities.”  http://icasualties.orj 


LED  Deaths  By  Year 


Hostile  Deaths 
m  IED 


2Qa7  ZQQp  2Qq<! 


Abbreviated  I  nj  ury  Scale 


TABLE  1.  Abbreviated  Injury  Scale  (A IS) 


AIS  score 

Injury 

1 

Minor 

2 

Moderate 

3 

Serious 

4 

Severe 

5 

Critical 

6 

Probably  lethal* 

*  Although  a  perfect  linear  correlation  with  an  AIS  ot  6  and  mortality  does 
not  exist,  survivability  is  unlikely. 


Examples  of  AIS  2 

■  Major  skin  laceration  or  avulsion  with  <20%  blood  loss 

■  Nerve  contusions  or  lacerations 

■  Vertebral  dislocation  without  fracture 

■  Herniated  disc  without  nerve  root  damage 
Lower  extremity  bone  fracture 


Center  for  Disease  Control  and  Prevention,  http://www.cdc.gov/mmwr/preview/mmwrhtml/figures/r801altl.gif,  accessed  on  April  30,  2010. 
Association  for  the  Advancement  of  Automotive  Medicine [1 990),  The  Abbreviated  Injury  Scale,  1990  Revision.  Des  Plaines,  IL. 
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Latin  Hypercube  Sampling 


Latin  Hypercube  Optimal  Latin  Hypercube 


http://people.sc.fsu.edu/-burkardt/m src/lcvt dataset/lcvt dataset.html,  accessed  on  December  5,  2009. 
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Model  Comparison 


Without  floorpad: 

500  data  points 

p  _  y  (  3.809— 0.03 95 4s£-j4  +  0.4289sc  +  0.003446npeaj^+0.0002 —  0.0 00001 

“neck  ^ 

F spine  =  383  -  462 sEA  +  416sc  +  lAapeak  +  262 sEAsc  +  0.7 sEAapeak  +  scapeak  -  232 sr2  -  0.0006a^afe 

Fcombined  tibia  ^7  3"  495sc  +  3 .7  Ct  peak  0.38sc(2pefj£  +  99sg  +  0.0003tfpg( 


Ri 

0.985 

0.979 

0.994 


With  floorpad: 

300  data  points 

u  —  7  (3.84  +  0.125^  +0.885c  +  0.002apeafe  +  0.058sEj45(:  +  0.0000845^ajjeafc-0.0000635capeafc-0.0585f>|-0.145c”0  00000054apeajt] 

**neck  —  e 

u  —  ?  (^.664 +  0. 12sea  +  0. 8 lsc  +  0.002 apeak  +  0.0625^ 5C  +  0.000087 sEAapeajc- -0.000068 scapeak- 0.0595^ -0.135^- 0.000000 56a^afc) 

** spine  ~  e 

Fcombined  tibia  =  332  -  245sc  -  80.23s/  +  1.3 apecLk  +  35.84scs/  +  14.0s/  +  0.0012aJeafe 


M 


R2 


0.952 

0.946 

0.976 


Optimal  Design  Laboratory 
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Formulation  1  Data 


Without  Floor  Foam  With  Floor  Foam 


Vehicle 
Mass  (kg) 

EA 

Stiffness 

Cushion 

Stiffness 

Floorpad 

Stiffness 

Maximum 
Injury  Ratio 

2000 

0.2500 

4.0000 

0.1000 

0.8219 

2500 

0.2500 

2.2457 

0.1000 

0.6125 

3000 

0.2500 

1 .6494 

0.7522 

0.4655 

3500 

0.2500 

1 .4056 

1.0632 

0.3783 

4000 

0.2500 

1.2530 

1 .2460 

0.3203 

4500 

0.2500 

1.1080 

1.3507 

0.2666 

5000 

0.2500 

1.0172 

1 .4267 

0.2351 

5500 

0.2500 

0.9392 

1 .4854 

0.2099 

6000 

0.2500 

0.8702 

1 .5704 

0.1892 

6500 

0.2500 

0.8120 

1 .8208 

0.1731 

7000 

0.2500 

0.7603 

1 .8867 

0.1598 

7500 

0.2500 

0.7166 

1 .9425 

0.1491 

8000 

0.2500 

0.6736 

1 .9975 

0.1393 

8500 

0.2500 

0.6410 

2.0391 

0.1323 

9000 

0.2500 

0.6042 

2.0860 

0.1247 

9500 

0.2500 

0.5777 

2.1199 

0.1195 

10000 

0.2500 

0.5449 

2.1619 

0.1134 

10500 

0.2500 

0.5286 

2.2053 

0.1104 

11000 

0.2500 

0.5036 

2.2145 

0.1060 

11500 

0.2500 

0.4804 

2.2441 

0.1022 

12000 

0.2500 

0.4588 

2.2716 

0.0986 

Vehicle 
Mass  (kg) 

EA 

Stiffness 

Cushion 

Stiffness 

Maximum 

Injury 

Ratio 

2000 

1.5000 

2.0000 

0.8616 

2500 

1.1082 

2.0000 

0.6597 

3000 

0.6323 

2.0000 

0.5175 

3500 

0.2962 

2.0000 

0.4328 

4000 

0.2500 

1.7909 

0.3757 

4500 

0.2500 

1.5907 

0.3333 

5000 

0.2500 

1.4406 

0.3029 

5500 

0.2500 

1.3039 

0.2760 

6000 

0.2500 

1.1905 

0.2543 

6500 

0.2500 

1.0942 

0.2362 

7000 

0.2500 

1.0110 

0.2208 

7500 

0.2500 

0.9415 

0.2081 

8000 

0.2500 

0.8789 

0.1968 

8500 

0.2504 

0.8183 

0.1860 

9000 

1.5000 

0.8094 

0.1720 

9500 

1.5000 

0.7929 

0.1641 

10000 

1.5000 

0.7754 

0.1558 

10500 

1.5000 

0.7582 

0.1476 

11000 

1.5000 

0.7437 

0.1408 

11500 

1.5000 

0.7284 

0.1336 

12000 

1.5000 

0.7178 

0.1286 
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Formulation  2  Data 


Without  Floor  Foam  With  Floor  Foam 


Vehicle 
Mass  (kg) 

EA 

Stiffness 

Cushion 

Stiffness 

Probability 
of  Failure 

2000 

1.5 

2.0 

4.60E-01 

2500 

1.5 

2.0 

2.45E-01 

3000 

1.5 

2.0 

9.93E-02 

3500 

1.5 

2.0 

2.97E-02 

4000 

1.5 

2.0 

6.43E-03 

4500 

1.5 

2.0 

9.90E-04 

5000 

1.5 

2.0 

1 .07E-04 

5500 

1.5 

2.0 

8.06E-06 

6000 

1.5 

2.0 

4.20E-07 

6500 

1.5 

2.0 

1.51E-08 

7000 

1.5 

2.0 

3.69E-10 

7500 

1.5 

2.0 

6.16E-12 

8000 

1.5 

2.0 

6.99E-14 

8500 

1.5 

2.0 

5.55E-16 

9000 

1.5 

2.0 

0.00E+00 

9500 

1.5 

2.0 

0.00E+00 

10000 

1.5 

2.0 

0.00E+00 

10500 

1.5 

2.0 

0.00E+00 

11000 

1.5 

2.0 

0.00E+00 

11500 

1.5 

2.0 

0.00E+00 

12000 

1.5 

2.0 

0.00E+00 

Vehicle 
Mass  (kg) 

EA 

Stiffness 

Cushion 

Stiffness 

Floorpad 

Stiffness 

Probability 
of  Failure 

2000 

1.65 

4.0 

0.10 

4.61  E-01 

2500 

1.65 

4.0 

0.10 

2.46E-01 

3000 

1.65 

4.0 

0.10 

1.00E-01 

3500 

1.65 

4.0 

0.10 

3.01  E-02 

4000 

1.65 

4.0 

0.10 

6.54E-03 

4500 

1.65 

4.0 

0.10 

1.01E-03 

5000 

1.65 

4.0 

0.10 

1.10E-04 

5500 

1.65 

4.0 

0.10 

8.36E-06 

6000 

1.65 

4.0 

0.10 

4.40E-07 

6500 

1.65 

4.0 

0.10 

1.59E-08 

7000 

1.65 

4.0 

0.10 

3.94E-10 

7500 

1.65 

4.0 

0.10 

6.66E-12 

8000 

1.65 

4.0 

0.10 

7.65E-14 

8500 

1.65 

4.0 

0.10 

5.55E-16 

9000 

1.65 

4.0 

0.10 

0.00E+00 

9500 

1.65 

4.0 

0.10 

0.00E+00 

10000 

1.65 

4.0 

0.10 

0.00E+00 

10500 

1.65 

4.0 

0.10 

0.00E+00 

11000 

1.65 

4.0 

0.10 

0.00E+00 

11500 

1.65 

4.0 

0.10 

0.00E+00 

12000 

1.65 

4.0 

0.10 

0.00E+00 

a 


'peak 


1756.7  G’s 


a 


peak 


1754.5  G’s 
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